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Summary 

Gallium nitride (GaN) is a most promising wide band-gap semiconductor for use in high-power 
microwave devices. It has functioned at 320 °C, and higher values are well within theoretical limits. By 
combining four devices, 20 W has been developed at X-band. GaN High Electron Mobility Transistors 
(HEMTs) are unique in that the two-dimensional electron gas (2DEG) is supported not by intentional 
doping, but instead by polarization charge developed at the interface between the bulk GaN region 
and the AlGaN epitaxial layer. The polarization charge is composed of two parts: spontaneous and 
piezoelectric. This behavior is unlike other semiconductors, and for that reason, no commercially 
available modeling software exists. The theme of this document is to develop a self-consistent approach 
to developing the pertinent equations to be solved. A Space Act Agreement, “Effects in AlGaN/GaN 
HEMT Semiconductors” with Silvaco Data Systems to implement this approach into their existing 
software for III-V semiconductors, is in place (summer of 2002). 


Introduction 

This report describes a method to model gallium nitride (GaN) with the Silvaco product ATLAS 
and its subroutine BLAZE. Unlike other semiconductor systems, GaN requires the addition of strains 
at all parts of the device to include the effects of piezoelectricity. One must also include a spontaneous 
polarization at certain interfaces. For these reasons this report consists of 1 1 sections. Section 1 introduces 
notations and definitions, in the area of linear elasticity. Section 2 discusses piezoelectricity and gives the 
complete set of basic equations to be solved in a general GaN device. Section 3 introduces spontaneous 
polarization. Section 4 indicates the additions needed in the Poisson-Schrodinger solver for GaN analysis. 

Section 5 gives some critical and basic details on the layer structure of AlGaN/GaN HEMTs. Section 6 
covers band-gap variations, while the remaining sections treat necessary issues to complete the analysis. 


Discussion 

Our goal is to model AlGaN/GaN HEMTs using Silvaco’s existing software platform. For the 
AlGaN/GaN system and other III-V nitrides, the effects of both spontaneous and piezoelectric 
polarization fields must be included in an analysis. These fields dominate the properties of the 2DEG, and 
this report serves as a starting point to incorporate them. The polarizations enter the problem via the 
constitutive equation relating D,E, and P : 


D = e 0 e r E + P 
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Where D is the electric flux density, E the electric field, and P the polarization of the material. 
This then influences Poisson’s equation as follows. Taking the divergence of the previous equation we 
have 


V • D =e Q e r V • £ + V • P 


Then Poisson’s equation is 


V-Z> = p 


or 


e o e r (y ' ( - V<j>))+ v P -P 
- e o e r V 2 ^p-V-P 

Thus the form using the electrostatic potential has an “extra” term due to polarizations not included 

in the relative dielectric constant e r . The added polarization P is composed of both the spontaneous and 
piezoelectric parts. The spontaneous part is modeled with an ab initio result of a theoretical calculation, 
while the piezoelectric portion must be determined from the stress/strain fields of the device. The 

constitutive equations for the stress/strain and D, E , and P are : 

S = d T -E + S e :T 
D = k s E + d:T=>P = d:T 


** ** r r 

The quantities T,S,S , d, k , are the stress, strain, compliance, piezoelectric, and permittivity 
tensors, respectively. 

The fact that P modifies Poisson’s equation means that it also modifies the energy levels for the 
2DEG as obtained from the Schrodinger equation. The strain also modifies the band gaps, via the 
deformation potential terms. 

This report consists of eleven sections. Section 1 introduces notations in the area of linear elasticity. 
Section 2 discusses piezoelectricity. Section 3 treats the concepts and properties of spontaneous 
polarization, as this is a dominant term, yet it is not well understood. Section 4 develops the forms for 
Poisson’s and Schrodinger’s equations. Here, the inclusion of dipoles and dangling bonds at the interface 
is addressed. Section 5 gives an account of the issues concerning the stresses and strains in AlGaN/GaN 
systems. Section 6 discusses the very important band-gap issues. Section 7 is just one page covering the 
variation in thermal conductivity with temperature. Section 8 gives numerical values for many of the 
terms discussed in other sections. It shows that the polarization terms are indeed dominant, as shown in 
published works. The effects of residual strains, strains due to interface mismatch, and oxide and metal 
pads are covered. Section 9 briefly visits some of the issues concerning heterojunctions. Section 10 gives 
data on many of the known traps discussed in the open literature. Finally, section 1 1 is a discussion of 
amplifiers and the best power levels reported to date. The appendix gives the closed-form expression for 
the field due to the spontaneous polarization. 
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1. Elasticity 


This section serves to introduce the needed definitions in elasticity theory that will be utilized later. 
Let u = mechanical displacement, then the strain is defined as (ref 1) 



f du , + du j 1 


l dr j dr > ) 


hj =x 9 y,z 


( 1 ) 


where S v means / = direction of displacement (or force) 

j = normal to surface (face) 

Strain is commonly given the symbols e^, w**, . . . 

Now 


S XX 

Sxy 

s x: ~ 


^1 



Sxy 

Syy 

S r- 

A 


*2 

2" ^4 

S x , 

Sr. 

S.. 


i^5 

4 

S 3 


( 2 ) 


Following Auld (ref 1), Nye (ref 2), and many others, we use the convention that compressive strain 
is negative, while tensile is positive. Note that the factors of Vi and sign convention vary in the literature. 
Write 


5 = 


Si 

s 2 

S 3 
*4 

S 6 


The stress is given as 


T 

XX 

T 

1 xy 

T 

1 xz 


'Ti 

t 6 

t 5 ~ 

T 

xy 

T 

J yy 

T 

- 

t 6 

Tz 

t 4 

T 

1 xz 

T 

J y= 

T zz 


Js 

t 4 



Alternate symbols for 7} are q^, . . . 


( 3 ) 


( 4 ) 
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Write 


T = 


Tj < 0 for compression 
Tj > 0 for tension 


( 5 ) 


Assume linear theory is applicable, then Hooke’s law is valid, so 


T ij ~ c ijkl S U 


( 6 ) 


or 

Stfkl =s ijkl^kl 

where c y */ are the stiffness constants (N/m 2 ) and 5 yW are the compliance constants (m 2 /N). 
The compact notation is 

T -c:S 
S = s:T 


( 7 ) 


(6a) 

(7a) 


where : stands for the tensor dyadic product. GaN and AIN form in either the Wurtzite (Wz)(hexagonal) 
or Zincblende (cubic) crystal classes. For Wz with symmetry notations 6 mm, Cg v , . . . 


To summarize 


Ti 

t 2 


c n 

c 12 

c 13 

0 

0 

0 

c 12 

c \\ 

c 13 

0 

0 

0 

c 13 

c 13 

c 33 

0 

0 

0 

0 

0 

0 

c 44 

0 

0 

0 

0 

0 

0 

c 44 

0 

0 

0 

0 

0 

0 

c \\ ~~ c \2 

2 

c u 

c 12 

c 13 

0 

0 

0 1 

c \2 

C 11 

c 13 

0 

0 

0 

c 13 

c 13 

c 33 

0 

0 

0 

0 

0 

0 

c 44 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

C U ~ C 12 


$5 

$6 


( 8 ) 


( 9 ) 
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*11 

*12 

*13 

0 

0 

0 

7 



* s 12 

*11 

*13 

0 

0 

0 

T 2 

s 3 


*13 

*13 

*33 

0 

0 

0 

7 

s 4 


0 

0 

0 

*44 

0 

0 

7 

*7 


0 

0 

0 

0 

*44 

0 

7 

.V 


0 

0 

0 

0 

0 

2(*n “* 12 ). 

.7. 


( 10 ) 


The determination of the five c v values is rather difficult. Presently GaN must be grown on a substrate 
(Sapphire, SiC, ZnO, Si, etc.). The epitaxial growth may be MOCVD, MBE, HVPE, among others; and 
the growth rates, gas pressures, and temperatures vary considerably. Under the present situation the films 
vary greatly. They differ in stoichiometry, defects, crystal faults, surface morphology, and film thickness 
(tun to pm). For these reasons, among others, one may only expect ranges of values for the fundamental 
constants for GaN or AIN. Table I gives two sets of values for the c v . The first column was derived from 
acoustic velocity measurements (ref. 3) on plasma-induced MBE films on Sapphire without a nucleation 
layer, with thickness between 800 and 1300 nm. The second column is from best guesses after averaging 
entries from the literature (refs. 4 to 8). 


TABLE I.— STIFFNESS CONSTANTS FOR GaN 


GaN 

c v in GPa 

c tJ (average, 
from literature) 

1 Difference 
(percent) | 


370 

382 

3.2 

C12 

145 

145 

0 

*13 

110 

106 

3.6 

*33 

390 

389 

0.26 

*44 

90 

103 

14 

*66 

112.5 

123 

9 


TABLE II.— STIFFNESS CONSTANTS FOR AIN 


AIN 

c tJ in GPa 

c tf (average, 
from literature) 

1 Difference 
(percent) | 

*11 

410 

390 

4.9 

*12 

140 

125 

10.7 

*13 

100 

120 

20 

*33 

390 

395 

1.28 

*44 

120 

118 

1.7 


Given the c tJ the s,j are (ref. 1 ) 


let 


C A ~ c 33 ( C 1 1 + c 12 ) “ 2 Cj 2 3 


then 


( 11 ) 




1 C 33 , 1 

2 |_ C A C 11 -c 12 _ 


( 12 ) 
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( 13 ) 



5 13 ~~ 


c 13 


(14) 


5 33 ~ 


C 11 +c 12 


(15) 


5 44 ~ 


1 

C 44 


(16) 


Table III gives the s v using best guess values for the c, r 


TABLE m.— COMPLIANCE COEFFICIENTS FOR GaN AND AIN. 



GaN 

AIN 

( 2 > 

s u ^UlO- 12 

V J 

3.175 

3.5319 

*12 

-1.044 

-1.0136 

*13 

-0.58072 

-0.76506 


2.887 

2.9965 (2.8 from experiment) 

*44 

9.7087 

8.474 


While GaN, AIN, and Ah Gai_* N are anisotropic, many papers treat them as isotropic for simplicity. For 
an isotropic material, we have 


X + 2p A, 

X A, + 2p 

A. X 

0 0 

0 0 

0 0 


X 0 0 0 

X 0 0 0 

X + 2p 0 0 0 

0 p 0 0 

0 0 p 0 

0 0 0 p 


(17) 


Where X and p are the Lame constants. Two other constants are normally used instead of X and p; 
namely E and v. Young’s modulus and Poisson’s ratio, respectively. They are given by 


E _ + 2p) 

(18) 

X + p 

X 

2{X + p) 

(19) 
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The bulk modulus is defined as 


The 77 S relationships in terms of E and v are (ref. 9) 

r “ = (UvTiTiv) [(> - v >*» + V (v + )] (21) 

T„ = (| + 4 - 2v) ~ ' V)S ” * + ' S " )1 <22 » 

T -, - ( Uv )(,_ 2v ) [<' - + v < 5 - + S » )1 (23) 


f 

T -T — C 

* ^(l+v) * 

(24) 

r 

T —T — 9 

*--~(l + v) ~ 

(25) 

T -T — ^ C 

** (l + v) v 

(26) 

S xx =^[T xx -v(T )y +T :: \ 

1 (27) 

S >y=\Vyy-'’{T xx+ T zz ) 

I (28) 

S==^K-* ( T * +T J 

(29) 

5 - 0 +v ) r 

(30) 

s _( 1 + v ) r 

£ *y 

(31) 

s = =(i+ / ,7 '“ 

(32) 
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There are several special cases: 
Plane strain (S :: = 0) 


Txx (l + vXl-2v)^ +V ‘ S »1 

| (33) 

Tyy ~ (l + vXl - 2v) ^ + V5xr 1 

(34) 

Ev r i 


^'-(l + v Xl-2v)^ + ^ 

(35) 

T - ^ $ 

*y i+ v v 

(36) 

yP 

II 

II 

O 

(37) 


s„=^p[(i-v)r„-vrj 

S.. = 0 


Sy: ~Sx: “0 

Plane stress (T z - = 0) 


T xx - - — 2 ( Sxx + vS yy ) 
1-v 


T - : 

yy j _ y 2 ryy 


(^>r + v $xx ) 


T.. =0 


T - — S 
*” 1 + v v 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
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T =T =0 

1 yz 1 xz w 


( 47 ) 


or 


- vT »’) 

S yy -^yy _v7 ^) 

= (l + v) 

*y ^ *y 

S,, = $* =0 

We may approximate the constants X, p, v, E, B, and M (the biaxial modulus) as follows. From 
equations (9) and (17) and table I, we see 

c, j = X + 2p = 370 = 382 
c 12 = A. = 145 = 145 

cj 3 = ^ = 110 = 106 

c 33 = ^ + 2p = 390 = 389 
C 44 = p = 90 = 103 
c 66 =p = 112.5 = 123 

From visual inspection it appears that a good choice is 

A. + 2p = 383 
and 

X = 126.5 

p = 128 


(48) 

(49) 

(50) 

(51) 

(52) 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
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With this choice, and using equations (1 8) to (20), we find (average values denoted with a carat ( A )) 


5 = 320 

(62) 

v = 0.25 

(63) 

5 = 212 

(64) 

E 

M = =427 

1-v 

(63a) 


In the literature I have found the following formulas 

B= = g_ (65) 

c n +c \2 + 2^33 - 4c 13 3(1 - 2v) 

E 9 

M = - = c n +c n 11 (66) 

1-v c 33 

Note, equations (65) and (66) permit the determination of 5 and v, as the c v are known, that is, two 
equations in the unknowns E and v. From the values in table I, we find 


Column 1 

5 =357 

v = 0.21 

5 = 206.6 

M = 453 

(67) 

Column 2 

5 = 377 

v = 0.20 

5 = 207.2 

M = 469 

(68) 


From the literature I have found the following reported values. 

5: 117, 160, 180, 188, 200, 204, 206,210, 220, 268 

These yield and average 5=195 GPa. 

v: 0.23, 0.25, 0.38, v = 0.29 

5:196,290 5 = 243 

M:413, 450, 479, M=447 
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From table I of reference 1 0 we fmd 


c n =363 
c n =128 

^13 =99 
c 33 = 3 8 5 
C 44 = 85 
5 = 195 
E = 273 
v = 0.27 


and the experimental values given (ref. 1 0 ) were 

B = 200 ± 20 
£ = 290 
v = 0.23 ±0.06 


( 73 ) 


(74) 


TABLE IV.— ELASTIC CONSTANTS FOR AIN FROM REF. 10 



From ref. 3 

Best guess 

I Difference 
(percent) | 

Cll 

410 

390 

4.9 

Cl 2 

140 

125 

10.7 

Cl3 

100 

120 

20 

C 33 

390 

395 

1.3 

C44 

120 

118 

1.7 


5 = 201 
v = 0.21 
M = 442 


The dispersion in the basic parameters reflects the wide variation in the films due to a variety of 
factors, as well as attempting to forge an approximate isotropic model. 

In most cases the films are assumed to be in “biaxial stress” which means 

Txx = Tyy, and T.. =0 (75) 


then equation (27) gives 


S xx E = T xx -vT }y =T xx -vT xx 


T — — S 

1 XX , XX 

1 -V 


(76) 
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which is apparently the reason for the definition of M in equation (66). The literature often uses the 
notation 


a xr “ -^OOOl e ;« (?7) 

where the subscript 0001 denotes the z-direction (parallel with the c-axis of the film). The crystal is then 
defined as having a Ga face. From equation (9) 

T 3 = ^3$] +C 13 S 2 +c 33^3 

If T 3 = 0, a “free surface,” then if Si = S 2 (biaxial strain), then 

0 = 2c 13 S' 1 + 


2ci3 

U- = u 


c 33 


(78) 


This result will be used often later. It can be shown that 


T = 

XX 


^ C \3 ~ C 33 ( C 1 1 + c 2 ) 


2 c 




(79) 


13 


where it is assumed that = 7^, and T :z = 0. This corresponds to equation (50) 

T„ 5.. 

2v " 

Using equations (60) and (61), X — 126.5, p = 128, then -El 2v = -639.2. Using the values in column 2 of 
table I, equation (79) becomes 


= -861 S z: (percent difference | = 26 percent 

This demonstrates the uncertainty in the magnitudes of the quantities in this section. 


(80) 


Estimation of Stress and Strain 

We assume the stress in the layers may be conceptually described by the superposition of three 
components: 

® TOTAL = ^ EXTERNAL + CT THERMAL + c INTRINSIC (81) 

The thermal part is the most straightforward, it is 

CT THERMAL =ct th =AaAT- (82) 

1 -V 
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Where Aa is the difference in thermal expansion coefficients (TEC), between the layers. Also AT is 
the difference between the growth temperature and the sample temperature (generally room temperature 
(RT)). The above assumes the substrate and epilayer are relaxed at the growth temperature (not true in 
general), and that no relaxation during cooling has taken place. 

The external portion (ctext) is due to loading by adjacent layers; lattice mismatch, metal pads, and 
any recognizable force. The last term is referred to as intrinsic or residual stress; it acts as a catch-all for 
stresses that do not fit into the two previous categories. The residual stresses are difficult to handle 
because they depend on growth parameters, which differ greatly for the basic growth methods. The basic 
methods are MBE, MOCVD, HVPE, and variations of them. Some of the parameters are growth 
temperature, pressure, gas flow rates, substrate crystal orientation, surface conditions, and preparation 
procedure. The GaN is often grown on a nucleation layer of either GaN or AIN. The properties of the 
nucleation layer include thickness, growth temperature, and doping. Finally, defects produced by 
vacancies, cracks, misfit dislocations, threading dislocations, etc. 

After the stress has been quantified the strain will still not be known exactly, as the material responds 
to stress in at least two ways. One is to strain via Hooke’s law, the other is to generate defects and 
propagate them throughout the material (relaxation). Finally, cracks do occur in GaN and sapphire, and 
many in the GaN are not visible by observation of the surface. 


External Stress/Strain 

Most often the strain due to lattice mismatch is assumed to be the dominant factor. Define the misfit 
parameter 


f m = 


a -On 


an 


(83) 


where ao is the relaxed lattice parameter and a is the strained value. The strain is then 

^xx~ —fm 

where the appropriate sign is chosen. Maeda (ref. 11) introduced a “lattice relaxation ratio, 7?” by 

*0 


(84) 


(85) 


This is very useful, as it permits some relaxation at the interfaces. The corresponding stresses are 

E Ac E 


2v c 0 


2v 




( 86 ) 


a xx 


E A a 
1-v a 0 


= M e. 


(87) 


where Ac/c 0 is the strain along the c-axis, and Aa/ao is that in the basal plane. The composite stress in a 
given film, or layers of film may be calculated at several levels of approximation. We start with Etzkom 
(ref. 12). 
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Figure 1 . — Epitaxial film of thickness hp on a substrate of thickness h s . Coordinate 
axes in the layers are z and t, respectively. Bending is due to differing thermal 
expansion coefficients. 


The figure is for the case when the net stress is compressive. The subscripts F an S represent film and 
substrate, respectively. The equations are as follows. 

Define the parameter 

r _ \Pf /( l~ v F )ty F _ M F h F 

[E s /(\-v s )]fi s M s h s 

(88) 

Assume the thermal strain dominates 


„ Act AT 
G F — M p 

1 + r 

(89) 

Act AT 

*s= Ms- r 

1 + r 

(90) 

^RT _ ( a F ~ a s) A ^ 
S 1 + r 

(91) 

_ RT ( a F~ a s) A T„ 
t c — r 

S 1 + r 

(92) 

h s -h F r 
2(1 + r) 

(93) 
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(94) 


R = 


J 

AaAT 


hi h F 1 \ 

6^ -2r) + 6(/7 F + 2r) + 2 5 + F ^ 


*F F)-f + ^ (95) 

Notice equation (95) gives the strain as a linear function of distance z from the interface. Equation 
(94) gives the radius of curvature. Observe that this formulation only requires material parameters and the 
thickness of the substrate and film. 

Next, the equation of Skromme (ref. 13): 


CT SURFACE ~ “ 


M s h\+M F h]r -3 Mphpjhs + *>f) 
6{h s + h F )h F R 


(96) 


where R is the radius of curvature which must be measured. Kozawa (ref. 14) has a slightly simpler form 
(assumes v the same in both layers) 


6(l — v)(/?£ + hp }hpR 


(97) 


TABLE V.— DATA FROM REF. 14 WITH OTHER 
ENTRIES IN PARENTHESES 



Sapphire 

GaN 

h F , h s (u) 

300 

2.5 to 48.8 

E (GPa) 

425 

196 (295) 

TEC ( 1 CT 6 ) 

7.5 

5.45 

V 

0.3 (0.25) 

0.3 (0.25) 


From Heame (ref. 15) we have Stoney’s formula: 

M s h] 

F 6 h F R 


(98) 


where for MOCVD at 1050 °C; sapphire: M = 450 or 479, depending on wafer; 6H-SiC: M = 586 GPa. 
Then 


o F =G xx =M F e xx (99) 

The composite stress was measured to be tensile with range 

0.14 < 0 ^ <0.29 GPa 


NASA/TM— 2003-2 1 1 983 


15 



The thermal component (by measurements) was -0.66 ± 1 GPa, and the theoretical thermal 
component is -1 .4 ± 0.7 GPa. Thus, the residual and perhaps some lattice portion dominates the stress 
in the given films. The radius of curvature R is determined by surface profilometry or laser reflection. 
Generally 


2 m < R < 20 m 

and in most cases R < 6 m. However, R = 0.58 m for HVPE films. 

Some experimental data (though incomplete) is helpful. From reference 1 6 a sample consisting of 
a 6H-SiC substrate with 3p of porous GaN, then 3p of HVPE GaN on top. The surface was in tension 
(wafer like a cup that holds water). The radius of curvature R was -6.5 m, the in-plane stress was 
0.64 GPa, and the strain along the c-axis was -1.12 x 10~ 3 . The substrate thickness was not given. 

A common range of parameters for GaN on sapphire is: 250 pm of sapphire with 500 A of AIN 
nucleation layer, and 1 to 5 pm of GaN. The residual strain in the GaN is of the order of 1 .0 x 10A 

General Case of Several Layers 

For the general case of several layers, one may use the results of Olsen and Ettenberg (ref. 1 7). The 
stress in any layer may be estimated from (see fig. 2) 




Figure 2. — General layer structure with thickness tj for 
each layer. Measured radius of curvature is R. 
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Here it is assumed that the elastic moduli (Young’s moduli) of all layers are the same and equal to E. In 
equation ( 100 )y is measured from the bottom of the 7 th layer. The terms in equation ( 100 ) are 


N 

'=!>, ( 101 ) 

;=1 


* = ' 3 6 II t itj e v 

(102) 

v ' =1 j> 1 y 


k,=1 i< j 

K ; =0 i = j 

K,. =— 1 i > j 

(103) 


In equation (102) the strain at the interface of the / -th and 7 th layer e v is related to the measured radius 
of curvature. However, it is assumed the strains e 0 are known, along with the thickness t, of each layer. 
If w~L , then one may assume spherical bending, then a two-dimensional (2D) solution may be obtained 
from the given one-dimensional (ID) by 


( \ CT yOOiD 

a '^D=^7 7 - 


(104) 


Reasonable choices for the interface strains are the larger of thermal, lattice mismatch, or residual 
strains (ref. 1 8). The reviews in references 1 9 to 22 should be consulted for more information. 


2. Piezoelectricity 

This section details equations we will use later. It is important to note that this discussion treats 
dielectrics, not semiconductors, that is, effects of carrier screening are ignored. We start with the defining 

equation for D, E, and P. 


D = e 0 E + P 


(105) 


Here 


P - P ELECTRONIC + PlONIC + POTHER ( 1 06) 

In most cases only the first two terms in equation (106) are assumed to exist, and Pother is assumed 
equal to zero. This case then has the constitutive assumption that the polarization is linearly related to the 
E field. 


Pelectronic + PionicA e 0 % e E 


(107) 
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where Xe is the polarizability. Then 


Z) = e 0 E+e 0 x e E 

= 6 0 e r & (^OTHER^O) 


which is the normal case. For anisotropic material we write 


~D X ' 


1 

* 

* 

* 

3 

X 

1 

'E x ~ 

Z)„ 

— 

K K,,. 

E , 

y 


yx yy y - 

y 

p. 


K :x K, y K„_ 

E 


We use k since e is used so often for strain. Thus 


— D = k-E 
e 0 


or 


— Dj = KyEj i,j = x,y,z 

e 0 

for GaN and AIN, we have (ref. 20) 

K n 0 0 

0 Kj j 0 
0 0 k 33 

for GaN, k u = 9.5 and k 33 = 1 0.4 
for AIN, Kj j = 9.0 and k 33 = 1 0.7 

Now when the dielectric material is also piezoelectric (following Auld, (ref. 1)) 

S = d T -E + s E \T 
D= g 0 k t -E + d :T 

which are the piezoelectric strain equations. Another form is 

T = -e:E + c E :S 
D — Gq k • E + € i S 


which are the corresponding stress equations. The superscripts mean constant electric field 
stress T , or constant strain S . By inspection we observe (now /other = ^p iez oelectric) 

P = d:f 


(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 

(113) 
, constant 

( 114 ) 
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where d v are the piezoelectric strain constants. We also have 


P = e:S 


( 115 ) 


where e v are the stress constants. For our Wz crystals we have (refs. 1 and 23) 




' 0 

0 

0 

0 

*15 

0" 

p y 

- 

0 

0 

0 

*15 

0 

0 

P. 


*31 

*31 

*33 

0 

0 

0 


Ti 

t 2 

h 

A 

t 5 

h 


(116) 


with the associations 7) = CTi = a*,, T 2 = a 2 = a >y , T 3 = g 3 = T 4 = <j 4 = g^, T 5 - ct 5 - g x: , T 6 - g 6 - a„, 
As in section 1, the range of the coefficients is large. 


The following strain constants are in units of 1 0 10 cm/V (see refs. 22 to 26): 


TABLE VI.— PIEZOELECTRIC STRAIN CONSTANTS OF GaN 



Measured 

Calculated 

Estimated 

Bulk 

Clamped 

Polycrystal 

d 3 , 

-0.9 

-1.16 

-1.7 

-1.9 

-1.4 

-1.0 

dx 


2.39 


3.7 

2.8 

2.0 

d ,5 

-3.1 

-2.75 






The following strain constants are in units of 1 0 10 cm/V : 


TABLE VII.— PIEZOELECTRIC STRAIN CONSTANTS OF AIN 



Measured 

Calculated 

Estimated 

Bulk 

Clamped 

Polycrystal 

<*3| 

-2.0 

-2.65 

-2.71 

-2.8 

-2.0 

-2.71 

d 33 

5.0 

5.53 

6.72 

5.6 

4.0 (3.2 film) 


^15 


- 3.4 




- 3.6 


Theory says d 3 \ ~ “(1/2) d 33 . The alternate representation is 

P = e:S 


or 




' 0 

0 

0 

0 

«15 

ol 

p y 

~ 

0 

0 

0 

e \5 

0 

0 

P. 


_«31 

e 31 

e 33 

0 

0 

°J 


( 117 ) 
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The following stress constants are in units of C/m 2 : 


TABLE vm.— PIEZOELECTRIC STRESS CONSTANTS OF GaN 



Literature 

Bulk 

Crystal 

Clamped 

Polycrystal 

eis 

-0.21 

-0.3 

-0.33 





^31 

-0.22 

-0.33 

-0.36 

-0.49 

-0.55 

-0.41 


-0.3 

e 33 

0.44 

0.65 

0.73 

1.0 

1.12 


0.85 

0.6 


The following stress constants are in units of C/m 2 : 


TABLE IX.— PIEZOELECTRIC STRESS CONSTANTS OF AIN 



Literature 

Bulk 

Crystal 

Clamped 


g 3l 


-0.6 



-0.5 


1.46 

1.5 


1.06 

0.85 



Also 


dj- • K — 


0 0 

0 0 

0 0 

0 d l5 

d l5 0 

0 0 



( 118 ) 


It can be shown that 


e Jj ~ c ji^Ij 


( 119 ) 
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Using this relationship we may determine the d,, 

^ 0 0 0 

<^31 =e 31 5 ll + e 32 5 21 +e 33 J 31 + e 34 ^41 + e 35 5 51 + e 36 5 61 

= e 3 j (sj 1 + s 21 )+ e 33 s 3 i = (- 0.35)(2. 131x1 0” 12 )+ (0.705)(- 5.8072 x 1 0 “ 12 ) 

= -1.16xlO~ 10 cm/V = -1.16xl0 -12 C/N 

^ 0 >j<0 ^ 0 

^33 =e 31‘ s 'l3 +e 32 5 23 +e 33‘ y 33 +e 34 5 43 +e 35 s 53 +e 36 s 63 

= 2e 31 5j3 +633533 = 2.39 xlO -12 C/N 

1 0 ^ 0 ^ 0 ^ 0 | 0 ^ 0 
^15 =e n s \5 +e 12 s 25 +e \3 s 35 + e \4 5 45 +e 15 s 55 +e \6 5 65 

= 6 15 5 44 = -2.75 xlO -12 C/N 

In our application, one often assumes the GaN to be completely relaxed, and the AlGaN to be in 
tensile biaxial strain. The resulting polarization field at the AlGaN/GaN interface is assumed to be as 
shown below. Start with 


P z = e 3\ u xx + e 3\ u yy> + e 33 u zz ( 120 ) 

Next assume a free surface (this implies the top of the AlGaN layer is free). This is a questionable 
assumption as the gate metallization loads the layer. Next assume biaxial strain, u ^ - u yy 


u zz 


c 33 


( 121 ) 


which is equation (78) of section 1. Then 


P-. =2e3\“ xx -2 — u xx e 33 = 2\ 

C 31 


e 3l ~ e 33 


£l3_ 

c 33 


( 122 ) 


The term in brackets has been measured to be —0.3 C/m 2 . An alternate representation starts with 


P : -^31°xx + ^3\ a yy +£ ^33 a :z 


(123) 


but 


^ XX ^1 1 ^ XX +Cy}Uyy +C\ 3 U., 


(124) 


and assume biaxial stress conditions hold ct av = a,,. Then 

P, =2 d 3x (c u +c n )u xx +d 3x c n u :: + d 33 cs : . (125) 

Next assume a free surface, so a zz = 0, 




-2 C ~^u x 

c 33 


(126) 
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then 


P : = 2d 3] (c n + c 12 >i xr +c, 3 


^-2— ^ 
V C 33 > 


' M x*^31 


(127) 


finally 


P, = 2d 


31 


C 11 + c 12 - 


4 

c 33 


(128) 


The corresponding electric field is (assuming D = 0, which means no charge is present) 

D _ 2J 31 


E. = — 


e n e 


0 t r 




"ll +c 12 _ 


4 


c 33 


(129) 


This field is about 6xl0 6 V/cm. equations 128 and 129 are found in the literature, (refs. 5, 27, and 28), 
and many have a factor of two before the c 13 / c 33 term. I interpret them as errors by attempting to add 

and Uyy, together as the “net strain” in the basal plane. For our case of Al* Ga^N, some authors 
multiply the P r in equation (128) by the mole fraction ;t to obtain the effective value. 

With the preceding concepts in place, we may introduce the entire set of equations to be solved in the 
GaN/AlGaN system. The first group consists of the mechanical equations of equilibrium (ref. 29) 


5cr S7Y,- QT 

w x t _j_ yx + kj± zx 


dx 


dr 


*y 


dy 

da 


dx 


+ - 


yy 




dy 

8 T, r 


d* n 

-qn— = 0 

ox 

zy 5<() 

— —-qn —-0 

dz dy 


dz 

dT. 


(130) 


yx 5 CT„ 0 <}> 

2 — + — —~qn— = 0 

dx dy dz dz 


+ - 


The last term in each equation is the “body force,” which we assume is the charge density times the 
electric field. Next, the electromechanical constitutive equations are 


°xx 


~c n C n C , 3 0 0 o' 

^ XX 


O 

o 

G yy 
a :: 


C 12 C n C 13 0 0 0 

C 13 C 13 C 33 0 0 0 

e_ 


0 0 e 31 

0 0 e 33 

d ^ J 
1 


0 0 0 C 44 0 0 

0 0 0 0 C 44 0 

0 0 0 0 0 — ^ 12 

yy- 

Jzx 


, 

o o o 
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° ° 

_ in _ 

° * ° 
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V u v U u 

2 J 
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'D x ~ 


’*11 

0 

0 " 

X’ 


' 0 

0 

0 

0 

0 

o' 

D y 

-G 0 

0 

*11 

0 

Ey 

+ 

0 

0 

0 

*15 

0 

0 

D. 


0 

0 

*33 _ 

E._ 


_*31 

*31 

*33 

0 

0 

0 


du _ dv 
dx ’ 8y 



=e„ 


dv dw du dw du dv 

dz dy ' - dz dx dy dx ^ 


( 132 ) 


(133) 


which couple the stress, strain, electric field, and electric displacement. The latter expressions are the 
compatibility expressions (ref. 29) relating strains to the material displacements. Notice the notation for 
mechanical displacements uses w, v and w ; rather than just u with a subscript. The last expressions to be 
used are the continuity, Gauss, and current equations. 



dJ x dJy dJ. n 
x + y + - =0 


a* 

dx dy dz 


<81 & 
1 

II 

J x = V-x n ( r ) E x+ D „ J ^ 

dx 

(134) 

b .- .a 

dz 

J y =V y n{r)E y +D„^ 

dz 


6D X 
~ + 

dx 

dD dD, 

, + ■ =An 

dy dz 

(135) 

n = n( 

r)=n(x,y,z) 

(136) 


cm 


where we have assumed just electrons are present. Similar terms for holes may be added for the most 
general case. Notice, we have 20 equations in 20 unknowns. This is an extension of the general elasticity 
problem to include a piezoelectric semiconductor. In elasticity theory, one has 1 5 equations in 
15 unknowns. 

The total set of equations has been provided for the sake of completeness. In principle they would 
be solved along with Poisson’s and Schrodinger’s equations in the Poisson-Schrodinger solver which 
presently exists in the software package. 
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3. Spontaneous Polarization 


There is some uncertainty about spontaneous polarization, Pspon , (see refs. 30 to 41). Bemardini 
et al., have performed ab initio calculations for spontaneous polarization in nitrides (refs. 30 to 34). 
Reference 34 discusses some of the controversial issues on the subject. On page 912 of (ref. 37) we quote 
“The latter distortion — henceforth called internal strain — determines a preferred polarity of the tetragonal 
axis and is responsible for the occurrence of spontaneous polarization.” I interpret the discussion to mean 
that there is an alteration of the position of atoms in a cell, which is deemed an internal strain (not 
apparent macroscopically), which causes the spontaneous polarization. From reference 38, “The so-called 
pyroelectric fields then arise from the temperature dependence of the spontaneous polarization when the 
temperature is changed.” Figure 1 in reference 39 shows both a macroscopic strain e and an internal 
strain w. Apparently u may exist when e = 0, so a polarization may exist due solely to u ; which would be 
a (macroscopically) strain-free polarization (spontaneous). A procedure to measure Pspon is given in 
reference 41. 

Pspon is screened by charges adsorbed onto the surfaces of the material, as well as being screened 
by free charge in the bulk. It is assumed (or defined) to be either uniform or zero in the bulk. One paper, 
not referenced here, stated it was zero in the bulk. Its presence is known only at interfaces with other 
materials (air being one of them). Pyroelectricity is the production of an electric field between two 
surfaces of the material held at different temperatures. 

P ? yro = P A T emperature (137) 

= [Px“x +P y a y + P=a : )&T (138) 


For GaN 


Px ~ Py 


= 10 4 


V 


m-K 


basal plane 


Pz 


7xl0 5 


V 

m • K 


along the c-axis 


Both Wz AljGaiJsf and GaN possess Pspon , and at their interface a positive bound surface charge exists 
(cubic GaN is piezoelectric but does not possess spontaneous polarization). This assumes a Ga face (only 
case we consider). This immobile charge is assumed (calculated) to reside in a layer 3 to 4 A thick, (I will 
draw it in the Al^Ga^N region for clarity). Figure 3 depicts two ideal samples with their bound 
polarization charges (no screening). 
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AIN 


(0001) 

C-axis 


p SPON 


+ + + + + + + + 

Pspon ( A,N * = -8X10“ 2 C/m 2 a z 


P SPON 

v 


+ + + 


GaN 


Pspon (GaN > = -2.9 xio - 2 c/m 2 a z 

Figure 3. — Ideal samples of AIN and GaN illustrating spontaneous polarizations. 


Note the z-axis and the c-axis are collinear. This relationship was assumed by Bemardini et al., and is 

important with respect to algebraic signs. Recall that the polarization vector P is defined to originate on 
negative bound charge and terminate on positive bound charge. The materials are assumed to grow from 
the bottom and increase in thickness as one moves up the c-axis (0001). From a simple charge super- 
position argument, a net positive bound charge density o(SURF ) exists at the surface. Analytically, we 
have 


c SURF =(8-2.9)xl° _2 C/m 2 = 5.1xl0~ 2 C/m 2 
= 3.19xl0 13 e/cm 2 

Recall that at a boundary (see fig. 4) 


Pn^Pn^-P^f 


® 


© 

Figure 4.— Definition for bound surface charge at 
interface between two polarized dielectrics. 

Here P 2 • n = P n2 , . . . 


For our case 

P„ 2 =P 2 -n= (- 8x 10 -2 a, )• a. =-8xlO" 2 C/m 2 
= ~P^-n =(-2.9xl0 _2 a.)-a. = -2.9x1 0" 2 C/m 2 
.-. P„ 2 -P^ =-8xl 0 -2 +2.9xl0 -2 = -5.1xl0 -2 C/m 2 
.•.-5.1xlO- 2 C/m 2 =-p^™ CE ) 

P £ 0 ^= a (SURF ) = 5.1x 10~ 2 C/m 2 


(139) 


(140) 
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Figure 5 summarizes the situation. 
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Al x Gai_ x N 


p SPON 
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©©©©© 

GaN 



Ga-Face 


Figure 5. — Ideal situation at AIN-GaN interface. Positive surface of AIN is reduced by negative charge 
on GaN surface. Result is a net positive charge at interface. 


Notice that the bound surface charge density is of the order of magnitude of the surface density of 
electrons in HEMT structures. 

Several papers, and in particular reference 40, have stated that the ab initio values should be reduced 
by at least a factor of 2. One reason put forth is that there may exist inversion domains at the interface, 
as well as islands of zincblende (cubic), observed experimentally, dispersed throughout the Wz phases. 
Other causes can of course be put forth; impurities and surface defects. It is also assumed (stated in one 
paper) that 


V-/>SPON=0 (141) 

as no bound volume charge density exists. This assumption seems reasonable if the internal strain is 
uniform and hence has no divergence. However, note that this assumption may be subject to future 
criticism. 


4. Poisson-Schrodinger Solver 

Figure 6 gives the coarse, initial band diagram in the region near the heterojunction. A key parameter 
is the step height A E c at the interface. Anderson’s model is presently used in some software, which gives 
the step the value 
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Al x Gai_ x N 


GaN 



2 


Figure 6. — Ideal triangular well for confinement of 2DEG at interface 
of Al x Ga«|_ x N and GaN layers. Conduction band edges are denoted 
as E c . 


“| XAl x Ga } _ x N ~'lGaN\ ( 142 ) 

where % is the electron affinity. We will discuss alternative ways to estimate A E c in later sections. The 
triangular well is composed of the step and the sloping side. The linear variation means a constant electric 
field is present. Apparently the establishment of the field comes about as follows. At the heterojunction, 
the step in E c along with the position of the Fermi level causes electrons to accumulate in the notch. The 
region into the bulk on the left is therefore depleted, and the existing donors establish the field. The field 
is modeled by the standard depletion approximation 

E(0) = l[N D =r i N s ] (143) 


where 

N d = depletion region sheet density 
Ns = free carrier sheet density 
r| = a weighting factor (~0.5 to 0.8) 

GaAs value at the heterojunction is 


2s(o)= 2.39x1 0 6 V/m (144) 

For our case of Al x Ga\_ x N/GaN, the situation is apparently quite different. From Sacconi (ref. 42), if 
we initially neglect all polarization, then the step height and Fermi levels (for typical doping) are not 
sufficient for a notch to form. Then electrons do not accumulate and create the 2DEG. After adding the 
spontaneous and piezoelectric factors, the notch forms, and the 2DEG is present. Apparently a depletion 
region does not form, as the polarization charge is sufficient to establish the field. If one does form, 
perhaps it is swamped by the polarization charge. One may look at the situation as follows: the 2DEG 
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forms to screen the bound spontaneous and piezoelectric charge residing at the interface. This introduces 
the question, “where do the electrons come from that establish the 2DEG?” There are valid reasons that 
state they cannot come from the bulk GaN, so apparently they come from donors somewhere in the 
Al x Ga^ x N region. This is not a trivial issue, as charge neutrality must hold, but the bookkeeping is not 
generally agreed upon in the literature. Instead of using the field expression developed from the depletion 
approximation, authors assume it should be represented by equation (129) in section 2, 


E. = 


Id 


31 




c n 

C U +c \2 

c 33 


(145) 


If the strain in the Al r Gai _ V N is tensile u xx > 0, then the field is oriented toward the bulk just as in the 
case of the depletion approximation. Figure 7 summarizes the discussion. 




Figure 7. — 2DEG formation, (a) By intentional doping, (b) By polarization charge. 


Poisson’s Equation 

Without spontaneous or piezoelectric polarization, Poisson’s equation is 


V 



N + D -N~ A +Nt -NJ +p-n] 


(146) 


where Np,N A , N, , N t are ionized donor/acceptor densities, and ionized donor/acceptor trap densities, 
respectively. Also 


e-e 0 e r 

<7 = |e| = +1.6 x10~ I9 C 


(147) 


The relative dielectric constant e r considers the standard electronic plus ionic polarizabilities. When 
both spontaneous and piezoelectric polarizations are included, the form is as follows. 

Start with 


V • Z) = p 
V-(s 0 E + p)=p 
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^ ' ( € O e r E + P piezo + E 5 PON ) = P 
V • [e (- V<t>)] + V • (/p [ EZ0 + /*spon ) = P 
€ [- V 2 <J)]— V(j)V 6 +V • (/> 3 IE 7 Q + ^SPON ) = P 

Neglect V e for simplicity 

V e— » 0 

V 2 <j) = — £ +— (v-/p ffizo +V-/spon) 

= - £ "-(p A (PIEZO) + p*(SPON)) 
e e 

then 

V 2 <t> = -^d-Na + K -N; + p-n\ 

E 

— “ [p* (PIEZO) + (SPON)] 
e 

We see that the potential now also depends on the divergence of the polarizations. Since 


(148) 

(149) 


(150) 


P =d:f = d :5 (151a) 

or 

P=e:S = e:u (151b) 

and the components of P are linearly related to the components of the stress or strain tensors; an accurate 
analytical model for the strain variation with position is mandatory. 


Schrodinger Equation 

Now we consider the effect of the polarization on the energy levels for the 2DEG. Assume the wave 
functions have the form 


M > J (x,y,z)=^(z)e^e jk ' x+jk ^ 


then write Schrodinger’s equation as 


n 2 6 


L*! 

l*, 2 +*f) 

2 dz 

v m* dz j 

1 2m* 


S, + V'{z%=Et, 


(152) 


(153) 
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We use an asterisk on the potential energy (PE) term V*(z ) , to denote that it represents the total potential 
energy field that establishes the well. The potential energy near the well is given in figure 8. The figure 
needs to be discussed in some detail. First, recall the relationships between the electric field E(z), the 
electrostatic potential <)>(z) and the electron potential energy band diagram (ref. 43) 



E c (GaN) 


z 1 

I 

I 

Figure 8. — First-order approximation of potential energy well at heteroj unction. 


E(z)=-^(z)=±^- (154) 

UZ Cf ctz 

Here E(z) is the magnitude of the component of electric field of interest, and E t the energy band edge. 

Here E t is a convenient reference energy level, and q is positive as given in equation (147). If only 
electrostatic forces are considered (the normal case in semiconductors), then the energy band edges take 
the shape of <| y(z) 9 with the negative sign treated appropriately. For heterojunctions, the band edges and 
§(z) do not track one another. This is due to the differences in electron affinities on either side, and how 
they change with position as one moves away from the interface. An electric field is generated with a 
magnitude of dyj(z)/dz. The variation of yjz) is difficult to determine thus some standard approximations 
have been used (ref. 44). Another force field is due to the spatial variation of the density of states (see 
refs. 45 to 47 for more information). The portion of the well near the notch is assumed to model the entire 
(or total) potential energy, and as one moves away, the field decreases and reaches flat-band deep into the 
GaN region. Most authors express V * (z) as 

V*(z) = -e$(z)+AE c U(z) (155) 

where U(z) is the unit step function. These two terms are sufficient to define the well. The first should 
include all charges in the vicinity, while the second is of a quantum nature. Another quantum term is the 
exchange-correlation potential V xc , but it is generally small, and may be neglected. Not all of the charge 
has been included in equations 146 or 1 50; namely, that which forms dangling bonds and interface 
dipoles. Some authors have also included the classic image-force potential. The exchange-correlation is 
(ref. 48) 


V xc =” 


1 + 0.7734* In 


HI 


V 7rar, 




(156) 
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where 


* 4n e 0 k/T 
a = - 


m* q 1 


K =s r (GaN)=e r (Al^Ga^N) 


r s= r s( z ) = 


n ( z ) 


-1-1/3 


(157) 

(158) 

(159) 


c(z) = - 


s ( z ) 


21 


(160) 


a = 


V 971 J 


(161) 


The classic image potential is 




8ti e 0 k a * 


v,m=- 


(k c -K b )q : 


(162) 


16;t e 0 k c (k c + k 6 )|z| 

where z is in the channel (GaN). The subscripts c and b stand for channel and barrier (AlGaN) regions. 

(k c -K b )q 2 


(163) 


V IM — 


1671 6 0 K c( K c +K *)(- Z ) 


(164) 


where z is in the barrier (z > 0 ). In practice V !M has been scaled by factors of 0.43, 0.53, and 0.75 for 
better fits to measured results. The dangling bonds have been represented by the ad hoc term (ref. 49) 


V db = 


47ie 0 e c d B 


( z < 0, channel) 


065) 


where e ( . = k ( . . The value for d K is a B J 3/4 where a B is the lattice constant in the barrier. In the barrier 


y db = 


4*e 0 e 6 d B 


(z > 0, barrier) 


(166) 
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where d c = a f V 3/4 where a c is the lattice constant in the channel (GaN region, z < 0). Apparently these 
are special cases of a more general form which varies with position. Assume 


V db {z) = 


"d<l 


4n& 0 e r (z 0 + \z\) 


(167) 


where n d is the number of dangling bonds at the interface. 

Since the junction is strained, the band gaps on either side will be perturbed via the deformation 
potential. Write the perturbation as (ref. 50) 


EE g ( z ) = “K | [«*r ( z ) + 11 yy i z ) + u -.-. ( z )] 


(168) 


where the negative sign reflects the experimental fact that hydrostatic compression ( u n < 0) causes the 
gap to widen. The deformation potential is a c . Therefore, from the existing literature, the most general 
expression for V* (z) is 

F*(z) = ^( 2 )+A£ g (z)+A£ f «(z)-e(|,(z) + ^(z) + F yw (r) + K rfB (z) (169) 

where E s (z) is the conduction band edge (ref. 50). This term is assumed piecewise constant on either side 
of the junction, and differs by the step height AE C . In effect, it reflects the difference in electron affinities 
between the two sides. The choice of assuming it is constant with position neglects the field due to d^dz 
(ref. 44). Only one paper included V db and many neglect V xc and V m as they are usually small with respect 
to A E c and <j >(z). Reference 50 is the only one to my knowledge to include E s (z) and A E g in V*(z). 

The previous paragraphs assumed no interface charge was present, and that <j)(z) is continuous across 
the interface. However, dipoles may exist at the junction (refs. 46, 51 to 58), and depending on the 
horizontal mesh size in a simulation, their effects should be included. Consider a 2D dipole formed by 
nearly adjacent layers of positive and negative charge as shown in figure 9. 


4>(z) 



A<Mz) < 




+ 

+ 

+ 

+ 

+ 



<Mz) 


d 


Figure 9. — Discontinuity in potential 
due to dipole (double layer). 
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The step change in <j>(r) is (ref. 59) 



e o 

where a = sheet charge density and d = separation between the sheets. 

If the horizontal step size in a simulation is greater than c/, then the dipole charge would be 
overlooked. Equation (170) allows its inclusion without an exceedingly small step (d is on the order of 
5 A). Some heterojunctions have negligible interface dipoles (ref. 58), but this may not be the case for 
the AlGaN/GaN system. Recall that a dipole (or double layer) is an assembly of charge as indicated 
above and d is less than other dimensions in the problem. Notice dipoles can add or subtract from A E c 
depending on their polarity. This may be used to adjust A E c for a given simulation. From experiments, 

A E c appears to vary from about 0.6 to 2.4 eV. Some of this variation may be due to variations in x(z), or 
doping, but a portion may indeed be due to dipoles. 

At the beginning of this section we showed the change in Poisson’s equation when polarization was 
included; now we show how it modifies V*(z). Rewrite equation (169): 

V*(z)=EM+^A z ) + ^Mzhe^{z)+yJz) + V IM (z) + V dh (z)+V DlP0LE (z)-eu p (z) (171) 

where we have added two more terms; V dipole and -eu p (z), the dipole and polarization potentials 
respectively. It is very important to notice that the -<t>„(z) now given the subscript n, in equation 1 71 is 
that given by equation 146, and not that given in equation 150. We do this to conceptually separate the 
effects of the normal charges included in Poisson’s right hand side, and the image, dangling bond, dipole, 
and polarization terms. 

The potential energy for a dipole is 


DIPOLE 


i z )~ “^DIPOLE 


where (ref. 60) 


^DIPOLE “ 


¥n 

4ti e 0 e r 


(172) 


(173) 


and vj/ = o T0T / dt 

a TOT ” charge on the positive side of the dipole 
di = separation between sheets 

Q = the solid angle subtended by the patch at the observation point (see fig. 10) 


Before we introduce polarization effects, recall the definition of the electrostatic potential 


Mr) 


1 f p(r')dv' 
4ne 0 |r-r'| 


(174) 
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Dipole (patch) 




Observation point 


Figure 10. — Geometry for determining potential of finite patch of double 
layer (dipole). 


When a polarization field exists, its potential is (ref. 61) 


471 Sf, * 


P(r')-(r -r']dv’ 


|F-Fl 3 


which may be written as 


U p 


00 


1 f P-Js' 1 f V'-iW 

4ne 0 ■k \r-r f \ 4^e 0 ^ v ' |r-r'| 

1 S Psb ds ’ , 1 f P vbdv' 

47i€ 0 ^|P — r'\ 47ie 0 *^ v r —r'\ 


(175) 


(176) 


where 


P sb= P ” = P n 

(178) 

Pv6=-V‘P 

(179) 


which are the bound surface and volume charge densities, respectively. The surface integrals are taken 
over surfaces bounding the polarization distribution. 

Now we attempt to ascertain the effects u p { r ) may have on the final energy levels for the 2DEG via 
V*(z). Apparently, the spontaneous polarization has the distribution of surface bound charge density as 
shown in figure 1 1 . 

The values of the charges are (from left to right) 

-8xl0 -2 —z-, 5xl0 -2 SL, and 2.9 xlO -2 ~ 
m 2 m 2 m 2 


These charges are and would contribute to u p ( r ) by the first term in equation 177. Assuming 
V • P s pon = 0? we have no further contributions from /spon • We neglect any pyroelectric terms, as no 
temperature difference can exist across 4 A. 

For the piezoelectric term, we have 

P ; d iezo = d : c = e:S (180) 

where the stress/strain may be decomposed into residual, thermal, and lattice mismatch components. 

It can also include stresses due to metal bonding pads. These terms will be dealt with in a later section. 
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For now, we can say that equation 171 gives the proper potential energy field and one must next 
evaluate u p { r ). 


Al x Gai _ x N/GaN 


GaN 


Substrate 



Figure 11 . — Assumed bound charge density at three main surfaces of Al x Ga 3 _ x N/GaN substrate HEMT. Charge 
is assumed to be developed by spontaneous polarization in Al x Ga-|_ x N and GaN layers. 


5. Material Configuration 

Figure 1 2 gives the general layer structure for our device. Our primary goal is to determine the strain 
at all parts of the structure. Starting at the nucleation layer we treat GaN/AlN/Al 2 0 3 . The AIN adheres to 
the sapphire and then allows the GaN to form on its surface. The AIN cracks and forms crystallites 
(sometimes) during its deposition. The GaN forms threading defects as it nucleates, and these threads 
spiral up through the material (never stopping) and continue to the AlGaN layer. It is unknown if new 
ones are initiated and propagate through the AlGaN. The defects have a density between 10 8 and 
10 12 cm -2 . The sheet value of 10 lo cm“ 2 corresponds to 2xlO l7 cnT 3 defects (assuming a sheet 5 A thick). 
Figure 13 gives a top view of the 2DEG and the negatively charged threads. The field of a charged thread 
varies from about 120 V/m to zero in a radius of 1 0 pm. Perhaps they should be modeled as filled traps; 
however, there is no obvious way to determine the equivalent trap density. 

Eastman’s research group at Cornell University state that a Ga face growth along the c-axis is 
guaranteed if an AIN buffer (nucleation) is used. The GaN is usually in compression and the A1 2 0 3 in 
tension. Figure 14 shows the stress at the interface. 

We neglect the fact that sapphire and SiC are piezoelectric, and only consider the strain in the GaN. 
First of all, the bound positive charge that supports the spontaneous polarization should exist at this 
interface. Most authors neglect this fact, and implicitly assume the charge is neutralized by immobile 
negative charge. This may be partially true, but some have noted a conducting channel at this interface. 
Eastman et al., observe “leakage” at this plane when using SiC and assume it starts to leak due to its 
relatively small band gap. In our modeling we may assume complete or partial neutralization with fixed 
charge with the remainder composed of mobile electrons. 

Often the strain at the GaN/sapphire interface is assumed to be due to lattice mismatch and/or TEC 
mismatch. Since the AIN deforms to relieve the ~16 percent lattice mismatch between GaN and sapphire, 
it is questionable if one can accurately calculate the stress or strain using the misfit parameter method. 
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Instead, we will use measured in-plane stress and strain. For reference, the misfit parameter for AIN on 
sapphire is (ref. 5) 


Gate 


Si 3 N 4 passivation 


Al x Gai_ x N 



* C-axis 
( 0001 ) 


z 

4 


Figure 12. — HEMT structure showing nearly all relevant portions. 


x 


2DEG 



Figure 13. — View looking toward GaN region from 
Al x Ga-|_ x N layer. 2DEG is assumed to be uniform 
with depleted regions around negatively charged 
threading defects. 



Figure 14. — Approximate stress profile near 
GaN/sapphire interface. 
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( 181 ) 


a 0 AIN = 3.1 12 A 
a 0 A1 2 0 3 =4.7589 A 


_ V3(3.112)-4.7589 
Jm ~ 3.112 
= 20.3 percent 


The VJ factor comes about because sapphire is rhombohedral, and the hexagons for the basal plane 
are rotated 30°. For GaN on AIN we find 


f m = 


3.1892-3.112 

3.1892 


2.4 percent 


^xx = -fm =-2.4 percent 


(182) 


and finally for GaN directly on sapphire we have (c-plane) 


or 


V3(3.1892)-4.7589 en , 

m 3.1892 


. 71(3.1892)- 4.7589 

f m = — 1 = 1 6 percent 

4.7589 


(183) 


This latter number is the one most often quoted, however, it is not really appropriate (ref. 5). These are 
the results from measurements: 


For GaN, the in-plane stress (refs. 62 to 65) = a n = a n 

0.2 GPa <|cr n | < 0.91 GPa 

From reference 15, |a n | = 0.66 ± 0.1 GPa. The lowest and highest values found in literature were 


I LOW 


= 0.18 GPa, ct, 


1 1 HIGH 


= 0.91 GPa 


The in-plane strain is w xv = u yy = u u =e n 


10 -4 <|m„|<5x10 -3 

The ranges found, including highest and lowest, are 

0.058 percents |m„|< 0.125 percent, | M ir | L0W = 0.016 percent, |w..| HIGH =0.57percent 
The calculated thermal stresses are GaN/sapphire Aa = 1.71 x 10^/K Ar~1030 K 

c TH = AaAT 300 = 0.7 GPa (184) 

lH 1-0.25 
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Recall equation 50 of section 1, and repeated after equation 79 reads 

T„=~S = (185) 

2v 

Using the measured a n and letting E = 300 GPa, we find 

e„= 0.33 percent 


which is within the experimental limits. 

For AljGai ,N, the measured stress and strains are 

0.8 GPa < <1.4 GPa 

0.01 6 GPa < Uja < 0.026 GPa 
0.0012 GPa < u,. <, 0.0057 GPa 

Moving up through the structure, we may state that the GaN is seldom intentionally doped, however, 
the minimum n-type carrier density is in the low 1 0 17 crrf 3 range. We have discussed the AEGa^N/GaN 
interface in section 4, and in general, the Al v Gai_ r N is in tensile strain. Most often, authors have assumed 
the GaN to be fully relaxed, which we will show is not true. The interface between the AUGa^N and the 
Si 3 N 4 passivation layer has not been reported on extensively, however, it is assumed to play an integral 
role in the problems presently plaguing the HEMT. It is assumed that the negative charge that supports 
the spontaneous polarization is neutralized by adsorbed ions. It is also assumed that donors may be 
present that actually deliver the charge that forms the 2DEG. There is disagreement as to the source of the 
2DEG charge. This interface is assumed by some to contain traps that interact with the gate, which causes 
the unexplained current collapse when radio frequency is applied. Some attribute the collapse phenomena 
to dielectric loss and bulk traps in the Si 3 N 4 . 


Stress/Strain With Position 

This is a most important topic, but unfortunately, not discussed much in the literature. Only one 
reference, Demarest (ref. 62), used the standard software package for stress/strain from ANSYS. The 
paper shows the strain variation around metal pads on Si02 layers. It should be mentioned that when GaN 
is grown on sapphire and SiC with AIN nucleation layers, polymorphs of Wz and zincblende regions may 
coexist. In some cases as much as 25 percent may be zincblende (on an Si wafer). During growth the GaN 
is in tension, but during cooling it goes into compression in the basal plane. Misfits form during the entire 
process, so the state of strain is highly variable from sample to sample due primarily to the different 
temperatures for growth, cooldown rates, and other factors. 

For reference, the basal plane strains at an interface are 

= ( 186 ) 

a SUfi 

where a S uB is the lattice constant in the layer in which the strain is sought. The corresponding strain along 
the c-axis is 
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For GaN, d and at, that is, as c decreases, a increases. The net strain is the combination of lattice 
mismatch, thermal mismatch (TEC), and residual strain. The latter often dominates, so we should rely on 
measured results as much as possible. For thin layers (up to about 2 pm), we may assume the stress 
is linear (ref. 63). Figure 15 depicts the measured stress for films of GaN on 6H-SiC (by MOCVD). In 
the region near 2 pm the material apparently relaxes and remains in horizontal tension for further growth. 
In figure 16 from (ref. 64), which shows the biaxial stress = a VJ , = a, , to be negative 

(compression). 



Figure 15. — Measured stress in GaN epilayer grown on 6H-SiC by MOCVD. Note, stress- 
free region occurs near 0.5 pm thickness (see ref. 66). 


£ 

<3 



t, gm 


20 


30 


Figure 16.— Reduction in stress with thickness of GaN on sapphire. 
Here = o-yy = an which is negative (film in compression). 
Absolute value is plotted. 


This trend is common; GaN on sapphire is often in compression while on SiC it is in tension. This, 
however, is not necessary, and the films may be in either state of strain on both substrates. It all depends 
on the mode, and temperature of growth, the nucleation layer preparation, the thickness of both the 
nucleation layer and the GaN layer, and the crystal orientation of the substrate and its surface condition. 

For GaN on sapphire (AI2O3), it appears that the strain decreases with thickness to some power, or 
almost linearly, or exponentially from the AIN layer (refs. 14, 65, 66). Table X gives the vertical strain 
with thickness t in pm (ref. 65). 
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TABLE X.— STRAIN ALONG THE c-AXIS 


VERSUS FILM THICKNESS 


f, pm 

Strain 

Fit equation result 

1 

12.3 

12.3 

5 

7.3 

7.3 

20 

4.6 

4.66 

60 

1.5 

3.26 


The fit equation is 


e„ = [l2.3 xlO -4 ]/ -0 32416 


(188) 


Table XI gives the basal plane compressive stress versus thickness (ref. 14), and it looks almost 
linear. 


TABLE XI.— MEASURED STRESS VERSUS 
FILM THICKNESS FROM REF. 14 


pm 

Compressive stress, GPA 

1 

-0.64 

10 

-0.59 

20 

-0.56 

50 

-0.47 


The paper by Detchprohm (ref. 67) shows the vertical (c-axis) strain varies exponentially from the 
nucleation surface. The strain does not relax until the layer is about 100 pm thick. This data was taken 
with MOC VD and HYPE grown samples (see fig. 1 7). 



Figure 17. — Relaxation of c-axis strain for GaN on sapphire. 
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The c-axis lattice constant was fitted to an exponential 


c(z)=5.185 + (0.006)e 00692 cinA,zinp 


( 189 ) 


then the strain with position is 


c(z)~ 5.185 
5.185 

0-006 ^-0.0692 
5.185 e 


Ac 

c _ 2^13 

Aa s n +s n 
a 


(190) 


(191) 


For Si 3 N 4 


-1.2GPa < a n <-0.02 GPa 


6. Band-Gap Issues 

The band diagram in figure 18 gives the ideal charge distributions at the Schottky barrier, the 
AlGaN/GaN interface, and the nucleation/substrate plane. 

We show positive charge on the Schottky barrier metallization that accumulates to neutralize the 
negative charge of the spontaneous polarization. The positive polarization charge and the 2DEG are 
shown at the heterojunction, and the spontaneous charge with some form of neutralization is shown at 
the nucleation layer. This is the minimum charge distribution for the device. Figure 1 9 shows the charge 
distributions for dipoles and dangling bonds that may exist at the Al^Ga^N/GaN interface. 



Substrate 




Figure 18. — Minimum charge distribution for an Al x Ga*|_ x N GaN HEMT. Charge at Schottky barrier, 
heterojunction, and GaN/substrate interface shown left to right. 
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Dipole that increases AE C 




Dipole that decreases AE C 



Dangling bonds 


Dangling bonds 


Figure 1 9. — Detail of possible charge 
conditions at heterojunction. 


The charge state in the Al^Ga^N and the Si 3 N 4 passivation layer are not considered here. The 
Al*Gai_,N barrier contains Al, but the total charge there is unknown. Often an Si donor layer is placed 
a few hundred angstroms from the interface to provide the 2DEG carriers. This doping apparently has 
a strong effect on the strain state of the barrier, which is a significant source of piezoelectric charge. 

We can omit the step discontinuity in <j> given in equation 1 70 of section 4, by combining it with the 

step; that is, eh<f> + A E c = A E' c . This keeps <(> continuous, if that causes problems with the software 
solvers. While this is convenient, it eliminates the force field of the dipole in the Schrodinger equation. 
Recall that the step AE C is due to the absence of available states in the barrier region, and not some form 
of repulsive field due to charge. The dangling bonds do not cause a jump in <)> and pose no problem. 

The range of A E c reported in the literature is about 

0. 1 7 eV < AE C < 1 ,8eV (192) 

Anderson’s model assumes 

= | XAljGa^N _ Xg3n| (193) 

However the uncertainties in x are very large (refs. 68 and 69) 

For GaN 2.1 <x <4.26 eV (194) 


For 


AljGaj.jN x(*) = 3-4;t for x< 0.75 
or = 3.25 -2.95 x x<0.75 


(195) 
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The steps in the band edges have been reported to be as follows 

A E v = 0.78x eV 

= 0.81 (x = l)and no strain 

v 7 (196) 

also 0.2 < AE V < 1 .4 eV 

For the nucleation sequence of GaN/AlN/Al 2 0 3 

A E v = 0.57 ± 0.22 eV for 5 to 20 A layers 

= 0.7 ± 0.24 eV (197) 

= 0.8 ±0.3 eV 


From these and knowing the band gaps, the step in E c is known. For the following range 

0. 1 5 < AE C < 0.56 eV 


the step is 


A E c ( x ) = 0.75[£ g (AljGa^N)- E g (GaN)J= 1 .96 jc 


For strain-free GaN, the following have been reported (see ref. 70) 


AE V = 1.36 ±0.07 


A E c 52 A E c 65 

— - = — or — - = — 
A E v 48 AE 20 


(198) 


(199) 


Thus, A E c is a fitting factor, and unfortunately it is one of the dominant terms in V*(z). Next we consider 
the change in E x (:) with strain. We write 

E g (STRAINED) = E g (BULK)± A E g (200) 

where E g (bulk) means the unstrained/undoped material. The change may be expressed as 

@ 01 ) 

where D i} is the deformation tensor and e tJ is the strain tensor. For GaN 


A = -0.88, A = 4.52, A = 5.4, A = -2.7, A = -2.85, A = 3.82 eV 


( 202 ) 


A simpler form using in-plane stress a„ = is 


AE 


g 


8Eq 

5a, i 


°ii 


( 203 ) 


where dE s /d(X\i = 21 to 27 meV/GPa 

= 42 meV/GPa for hydrostatic stress (and Si doped) 
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Sometimes one finds (ref. 71) 


dE 


1 = 42 meV/GPa 


d Pressure 

Other forms using both in-plane and c-axis strains are as follows: 


E (STRAINED) = E (RELAXED)h 


dE n 


de„ 




where 


dE c 33 

de :: 4 c 13 



GaN 


AIN 

D 3 (eV) 

5.8 


8.84 

£>4 (eV) 

-3.25 


-3.92 

2) is 





A^HYDRO 

+ A^SHEAR 

A^HYDRO 

= ~ a [ e xx 

+ G xr 

+2 e z: 

A^SHEAR 

= -b [^xx 

+ e xr 

-2 


(204) 


(205) 


(206) 


(207a) 

(207b) 

(208) 


where a = -8.16 eV, and b = 3.71 eV. In some cases the hydrostatic and uniaxial strains are given by 


= HYDRO" 


e 11 +- 


= UNI AX 


~T( e X _e ll) 


(209) 


along [OOOl] 


Another form is 


A E g =C e„ 


( 210 ) 


where C was determined by experiment to be 12 eV. 
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Finally, one finds in terms of stress 


de„ 


= 10.1-15.4 eV 


in terms of strain 


— = -6.1 eV e u < 0.0032 

11 

= -8.2 to - 13 eV (TENSILE) 
= -16 eV e n > 0.0032 


( 211 ) 


( 212 ) 


The expression for A E g in equation 168 in section 4 is that due to hydrostatic stress only. 

In figure 20 (based on a figure in the literature); we show the energy gap of GaN with in-plane strain. 
Some authors claim that the curve is symmetrical, that is, the gap increases for both compression and 
tension. 


4.0 


$ 3.8 

d 

<0 

o> 

§ 

<D 

iS 3.6 


’ 0.00 0.01 0.02 0.03 0.04 0.05 

€ n 

Figure 20.— Variation of bandgap of GaN with basal plane strain eu = €** = Cyy 
(see ref. 73). 


Now we consider the change in band gaps with temperature (refs. 73 to 75). 
For GaN 



ji2 

E e = 3.503 + 5.08 x 10" 4 eV 

g T - 996 K 

T 2 

or E„ =3.509 -7.32x1 0" 4 eV 

8 T + 700 K 

y 2 

or E e =£ p (r = 0)-8.32xl0 -4 — — eV 

8 gK 7 T + 836 K 


(213) 
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For AIN 


7-2 

E a = 6.118-1. 8xl0~ 3 eV 

g T + 1462K 

7-2 

E p (x, T ) = (3.42jc + 0.068)- 8.73 x 1(T 4 eV 

' v 7 T + 830K 

The band gap of Al^Ga^N is assumed to vaiy with mole-fraction jc as 

E g (x)=xE g (AIN) + (l - x)E g (GaN) - bx{\ - x) 

where the bowing parameter “6” varies from 0.53 to 2.6 eV in the literature. Some have even stated it is 0, 
±0.2, or ±0.3 eV. The following forms are found. 

E„- 3.43 + 2, + 0.7,> (2I6) 

= 3.43 + 2.56* 


(214) 


(215) 


7. Thermal Conductivity 


Let k(T 0 ) be the thermal conductivity at T 0 (normally 20 °C). Then the variation with temperature is 


k{t) = k(t 0 



W 

cm K 


(217) 


TABLE XB.— SUMMARIZES THE VALUES 
FOR OUR MATERIALS 


Material 

K 

r 

Sapphire 

0.25 to 0.50 

0.544 

SiC 

3.3 SI, 0.5 

0.524 

GaN 

1.3 to 1.7 


Si 3 N 4 

0.37 


AIN 

2.0 


Si 

2.0 



Some experimental results for AIN are below. 


TABLE xm.— THERMAL CONDUCTIVITY 
VARIATION WITH TEMPERATURE FOR AIN 


T 

Experimental value for k 

300 

2.85 

400 

1.80 

600 

0.96 

1000 

0.48 
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8. Orders of Magnitude of Terms and Approximations 


The exact form for the potential of a rectangle of charge is given in the appendix. It is tedious so we 
will assume the potential of a disc of charge is a reasonable approximation. The potential of a disc of 
radius “a” is 


<t>(z) = 




-V^7?-z 2 ‘ 


(218) 


where a 0 is the surface charge density. Assume the source-drain spacing is 2 pm, then a - 1 pm. For the 
spontaneous polarization at the interface, the effective charge density is 




= 5.1xl(T 2 



Evaluate at 70 A from the surface 


<(>(70A)=320V 


(219) 


( 220 ) 


This would be the value if no screening existed due to mobile charge. We will assume the screening can 
be modeled using the classic Debye-screening factor 

Ids ~ (221) 


where £& is the Debye length 


iD _ [Ik 

» 2p ° q ( 222 ) 

Q 

Po =n 0 q = \.6x\0 6 — r 
m 

The densities in the 2DEG are on the order of 10 19 cm -3 , which yields a length of about 8.25 A. The 
evaluating point of 70 A is then about 8.5 Debye lengths, and the reduction factor becomes 

f DS = 2x10^ (223) 

Thus, the potential is 64 mV, which is well within acceptable limits. 

Now we estimate the terms in equation 171 in section 4, to ascertain their relative magnitudes. The 
step in the conduction band is bounded as 


0.17eV <AE C <1.8eV (224) 

Therefore, all the other terms should be less than an eV or so. The change in the bandgap due to strain, 
A£ g , has a large dispersion. From section 6, we find 

20 meV < A E g < 230 meV (225) 


NAS ATM— 2003-2 1 1 983 


47 


While large, it is definitely within acceptable limits. For the electrostatic potential, we start with the 
electric field from equation 129 of section 2. Its value is 4.8x1 0 8 V/m, which is about 100 times larger 
than that for GaAs. It is very near the breakdown value (5 to 6x 1 0 8 V/m). At 70 A we have 

<J>(70A) = E s (o) (70A) = (4.8 x 1 0 8 V/m)(70xl0" 10 m) 

= 3.36 V .\<?<|> = -3.36eV 

which is high. If we reduce the field by a factor of 10, which is reasonable, then the value is acceptable. 
There are two justifications for the reduction by the factor of 10; one is the uncertainty in the stiffness 
parameters while the other is the chosen order of magnitude of the strain. The exchange correlation is 
-38.6 meV and the image force gives -0. 146 meV. The dangling bond gives 

V db = 20 Jn d meV (227) 

where n d is the number of bonds. This would give a very large value for even one bond. Using the 
screening factor reduces it by 10^, which brings it into a reasonable range. The dipole contribution gives 
-896 eV without screening (for just 100 electrons on a sheet). The screening reduces it to -89.6 meV, 
which is excellent. 

Now we estimate some of the strains that will contribute to Poisson’s equation 150 in section 4. 

We start with the strains due to metal pads and dielectric films on the Al*Gai_*N top surface. Consider a 
semi-infinite film of thickness h on a substrate. The stresses due to this “single edge” (refs. 76 and 77) 

(see fig. 21). 


T l = CT x* =- 


-2 h 


( 


Vi ) 


(228) 



Figure 21 . — Geometry for calculation of stress in region z > 0. Single 
stage at x = 0 develops stresses for x > 0. 
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T 2 =°yy = 


_ ~ 2v/ £ 


7 i — a — — 


f .. ^ 


v r i y 


-2A 




V n y 


r 4 = a v _ =0 


^5 - CT .tr - 


- 2/a' 


^ y2 ^ 

X, Z 


V '1 / 


T 6 = °.vv = 0 


(229) 

(230) 

(231) 

(232) 

(233) 


These are correct for x>0. They are not correct for large negative x (which is out of our range). 
When two edges are used we have the following situation (see fig. 22). 



Figure 22. — Geometry for loading due to pad of width Lg. 
Modeled as superposition of two single edges at x = 0 
and x = Lg. 


f E =°/dj+° g d g 


where a j,ts g are stresses produced at z = 0. 


T\ =°xx 


Ve 




(234) 


(235) 
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(236) 

(237) 

(238) 

(239) 

(240) 

(241) 

(242) 

(243) 

(244) 

(245) 

(246) 






( 247 ) 


Notice the coordinate system is not that used in other sections. To conform to that system observe the 
following changes: z (here) goes to (-z) there, and>> (here) goes to (-^) there. Thus, the only changes 
required are the signs of y andz. This has already been done in equations 245 and 247 above. 


NASA/TM — 2003-2 1 1983 


50 


A previous study using the above equations has been done in the AUGai^As/GaAs heterojunction 
system (ref. 78). It was observed that the stress changed from tension to compression at the mole-fraction 
x = 0.03. The crystal type was zincblende (cubic) so 


u 


n 


1 

j 


'o 

0 

0 

*14 

0 

Pi 

= 

0 

0 

0 

0 

e 14 

A. 


0 

0 

0 

0 

0 


0 

0 

e l4 


u 22 
w 33 
u 44 

u 55 


\_ u 66 


(248) 


Note, only shearing stresses cause P in this system. If grown in the [100] direction, then P = 0 ; if grown 
in any other direction, then P * 0 . For GaAs 

v = 0.23 

d u =2.6x1 0 -10 ^- (close to GaN) (249) 

P x = d u a, x 


P y = d X4 o yz 


(250) 


P. 



(251) 


Figure 23 shows qualitatively the differences between the “edge” model and a more accurate 
calculation using ANSYS (ref 79). Figure 24 is from (ref 80), which shows the piezoelectric charge near 
a gate pad. Observe the charge density oscillates from positive to negative as one traverses the device on a 
horizontal line. The densities are as large as 10 17 cm -3 , which is the same order of magnitude as the 
background density obtained in GaN. 
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Charge density, charge/cm 3 



Figure 23.— Doping profile and piezoelectric charge density p p2 calculated by "edge" load 
line method, and more accurate ANSYS finite element method (FEM) solution. 
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(a) (b) 



Figure 24. — Characteristic "bifoleum loops" of piezoelectric charge generated by edges of pad on surface. 
Numbers represent charge density; 17 means 10 17 crrr 3 . (a) [01 1]-oriented field-effect transisters (FET) on 
[100] plane : (b) [Oil]-oriented FET on [011] plane, (c) [21 1] -oriented FET on [ill] Ga plane, (d) [011]-oriented 
FET on [1 1 1] Ga plane. 


9. Heterojunctions 

In other sections various properties about heterojunctions were mentioned. This section seeks to 
gather much of the information I have gathered that may prove useful. First of all, the assumption that the 
interface strain is equal to that of the misfit parameter needs to be discussed. While it may serve as a 
guide, its reliability is questionable. First of all, the intrinsic and thermal strain components must also be 
calculated to estimate the strain at a junction. The most common interface is that of GaN on sapphire. The 
misfit parameters in GaN and sapphire are 


For GaN 


fm 


3.189V3-4.7589 

3.189 


= 24 percent 


For sapphire 


fm = 


3. 189V3 -4.7589 .. 

= 1 6 percent 

4.7589 


(252) 
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It would appear that the GaN would be in lateral tension (and thus in vertical compression), and 
visa-versa for the sapphire. Experimentally, while GaN is growing on sapphire, it is in tension. The 
corresponding TEC are 

For GaN a = 5.59x1 (T 6 /k 

For sapphire a = 7.3 x 1 0 6 /K 

This implies the GaN would be in compression. Theory would say the thermal stress is 

i | A AT ,£(GaN) 

Pu = AaAT — 1 - 

1 - v 

= (l .71 x 1 0~ 6 )(l 020) ^ 3 q 25 GPa (253) 

= 0.74 GPa 

Experimentally, the net stress is near 0.66 GPa, compression. 

The next combination is AIN on sapphire (which is the standard nucleation layer). The misfit 
parameters are 


3.189-4.7589 

For AIN /_ = = 53 percent 

m 3.112 


3.112-4.7589 

For sapphire f m = = 35 percent (254) 


The TEC for sapphire is a = 7.3x1 0 -6 /°C. It would appear that the AIN is in lateral tension latticewise and 
compression from TECs. 

Next we consider GaN/AIN 


. 3.189-3.112 ^ A . 3.189-3.112 

/_, - = 2.4 percent, /_ = = 2.5 percent 

3.189 m 3.112 


(255) 


The GaN should be in lateral compression latticewise but in tension from the TECs. When joined, the 
measured lattice is 3.1 84 A. Experimentally, the GaN is in compression. 

For GaN/6H-SiC 


fm = 


3.189-3.08 


= 3.4 percent 


3.189 

a(6H -SiC) = 4.2x1 0 _6 /°C 


(256) 


The GaN in compression latticewise and tension TEC-wise. In many experiments it is in tension, from 
0.1 to 0.2 GPA (for 3 -pm films). 
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For GaN/4H-SiC 


fm =' 


3.189-4.23 


= 33 percent 


3.189 

a(4H-SiC)=4.2xlO“ 6 /°C 


Both imply tension for GaN. 

The case for AfGa^N/GaN, the lattice constant for AfGa^N is 

A(x) = (-.0 77x + 3.189)A 


(257) 


(258) 


Then its misfit parameter on GaN is 


fm = 


-0.077a: + 3. 189 -3. 189 
-0.077a: + 3. 189 


(259) 


The TEC for AIN is about 5.3, as compared to GaN’s 5.59, so assume no added strains during 
cooling. Experimentally, AfGai^N is in tension. The in-plane stress is between 0.82 and 1.33 GPa. The 
in-plane strain ejj is between 0.016 and 0.026. Table XIII gives the vertical strain with mole-fraction 


TABLE Xin.— THE VERTICAL STRAIN WITH MOLE-FRACTION IN AfGa^N 


X 

0.05 

0.1 

0.15 

0.20 

0.25 

Gzz (percent) 

| 0.12 

0.24 

0.35 

0.48 

0.57 


The c-axis lattice constant in GaN is 

C(x)=(- 0.2661 x + 5. 17)A 0<a:<0.9 
The total spontaneous polarization is 

^spon(*)= -0.029(1 -*)- 0.08 lx 

m 

= (- 0.029- 0.052*) => 3.25* (lO 1 3 ) electrons/cm 2 

/>(PIEZO)GaN = 0.0163 a:—- 

m 

P (total) GaN = P s P0N (GaN ) + P piezo (GaN) 

= (-0.029 + 0.0163a:) 

nr 


(260) 


( 261 ) 


(262) 


(263) 
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^piezo (Al_ r Ga 1 _ ;t N)-4.26x(lO 2 ) — — 

m 

= (2.66 - 3X*)(l 0 1 3 ) electrons/cm ‘ 


= 2 




L33 


(*) 


P (totalXboth materials) = (- 8.4x - 1.9x 2 )x 1 0 6 


C 


= -0.096x 


m 


2.6 xl0~ 2 < /’(total) < 4 x 10 -2 

m 2 

The stiffness coefficients are 

c 13 (x)=5x + 103 = llx + 105GPa 
c 33 (x) = — 32x + 405 GPa 
The stress coefficients are (refs. 5, 81, 82) 


e 3 1 (x) = _ 0. 1 lx - 0.49 ■ 


m 


e 33 (*) = -0.73x + 0.73 —j 


m 


The relative dielectric constant is 

k(x) = -0.5x + 9.5 

The Schottky barrier heights are (ref. 83) 

e$ b =(l.3x + 0.84)eV Ni 
e<t> 6 = (2.44x + 0.9l)eV Ti/Al 

The effective masses are 

m e (x) = [o.3x + 0.2 (l - x)]w 0 
m hh (*) = [3.53x + 1 .76(l - x)]/n 0 

The electric field is 


m 


(264) 

(264a) 

(265) 

(266) 

(267) 

(268) 

(269) 

(270) 

(271) 

(272) 

(273) 

(274) 


NASA/TM— 2003-2 1 1983 


56 


(275a) 


t(x) = (7.9 x 1 0 -4 + 1 . 1 5 x 10 7 x)v/cm 


ore(;c)=(9.5;c + 2.1;t 2 )xl0 6 V/cm (275b) 

or 8(x)=1.2a:x 10 6 V/cm (275c) 

Finally, the surface sheet density is 

CTtot(*)= 0 - 068 * -% (276) 

rrr 

/7 5 (x) = 6a:x 10 12 — (277a) 
cm 

or w 5 (jc)= 5.5xx 10 13 --- - - (277b) 

cm 2 


10. Defects 

For unintentionally doped AlGaN the carrier density is near 1.0 17 cm 3 . The carriers are assumed to be 
caused by N vacancies, but there is no general agreement. Listed are some speculations on defects. 

1 . Surface traps cause the radio frequency compression. These are emptied and filled by the gate 
voltage. Some speculate that the gate voltage alters the strain in the barrier and passivation 
layers and thus the amount of trapped charge. 

2. Shallow traps cause I-V dc current loops. 

3. Deep traps in the buffer are the cause for the I-V dc current collapse. Hot electrons are 
introduced into the buffer for drain voltages above 20 V. Gate voltages have no effect on this 
process. 

A surface donor at 1 .4 eV, 1 .65 eV. 

A deep trap at 1.91 eV. 

A shallow trap at .95 to 1.0 eV with capture cross-section 4xl0 _12 cm 2 . 

For mole-fractions, the density of these traps are: 


TABLE XIV. — APPROXIMATE TRAP 


DENSITY IN AfGa^N 

F XJ 


X 

N, cm 3 

0.05 

5x 10* ' 

0.15 

2xl0 111 

0.25 

5x10 11! 


For semi-insulating (SI) GaN a trap exists at 2.15 eV below E c with density between 3x 10 9 and 
3xlO lo cnT 2 . For normal GaN, traps at 1 .8 and 2.85 eV below E c . A deep acceptor trap exists at 2.2 eV 
below E c with density ~1 x 10 10 cm -2 . 
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The nitrogen vacancy has a state 40 meV below E c , and the £> trap exists 0.55 eV below E c . The trap 
that causes dc current collapse is 1.19 eV above Ey. The dominant Mg trap is 150 to 200 meV above E v . 
A trap exists at 0.95 eV above E v and the characteristic yellow emission is due to a transition between 
levels 1 eV below E c to one 0.2 eV above E v . The following list gives the levels (I assume with respect 
to E v ) of many of the standard impurities found in GaN (ref. 22). 

Zn 0.34, 1.8-1 .9, 2.2, 2.5-2.6, 2.8 

Mg 0.25, 2.95, 3.2 (note, 0.15 to 0.2 in the paragraph above) 

Cd 0.55, 2.7, 2.85 

Be 0.7, 2.2 

Hg 0.41,2.9,2.43 

C 0.86,2.15 

Li 0.75, 2.23 

P 2.85 

As 2.58 

Si 0.225 


11. Amplifiers 

The best ever performance: 9.1 W/mm at 10 GHz on 4H-SiC (refs. 84 to 86). The maximum available 
is 12 W/mm, for then the channel goes above 300 °C and/or the SiC substrate starts to leak, and large 
anomalous current starts. 

Devices have been plagued by both dc and tuning-induced failures. 

Eastman says degrade most theoretical f by 30 percent. SiC has maximum thermal dissipation of 
27.8 W/mm, but starts leaking at 10.9 W/mm. SiC leaks more than sapphire due to smaller bandgap. 

Eastman 0.2 < x < 0.4 

Channel temperature < 300 °C 

P(max) at 10 GHz = 12.5 W/mm 

The channel temperature is 


^CHANNEL - To (SUBSTRATE) + R lh I D V DS 

R ,h = 24 [1 + (r - r 0 )x io _J 

Wnel(MAX)=500K 

The devices are predominantly thermally limited, so thermal simulations would be important. 
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Appendix 


z 



Figure 25. — Epitaxial Al x Ga-|_ x N layer (z > 0) on GaN (z < 0) with spontaneous 
polarization P and bound charge on plane z = 0 and z = z 0 . 


The field due to the spontaneous polarization of the epitaxial layer of Al^Ga^N is developed. 
Figure 25 gives the geometry. The spontaneous polarization is specified by 


P SPON ~ ~P ) a : (A1 ) 

Where Pq = 0.08 C/m 2 on the top surface and 0.051 C/m 2 on the bottom surface. This assumes no 
neutralization of the negative charge on the top surface due to adsorbed charge. It also models the net 
charge at the heterojunction. Starting with the first term of equation 1 76 of section 4 


u P (r)= <f 


P- fids' 


4tc < 


r-r' 


(A2) 


On the upper surface ft = a . , so P ■ n = -P 0 . While on the lower surface we have n - -a z , P ■ n = P 0 

Thus, the two surfaces at z = 0 and z = z 0 are just uniform charge sheets. 

On the lower surface we have 


r' = x'a x +y'a y 

(A3) 

r = xa x + ya v + za. 

(A4) 

I f - r'\ = yl(x - xf +(y-y'Y + z 2 ds’ = dx'dy’ 

(A5) 
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Using elementary steps and reference 87 we find 


r*o rvo 1 

P 0 dy'dx' 

° Jo 471 € 0 l(x-x') 

U(y + y'f + z 2 ] n 


From no. 200.01 (reference 87), we have 

P 


p ' An € 

Pq 
An e 


( r )= f°4 ^ _ ^o) + ^(y-yo) 2 +A ~~ 

i, 


dx' 


p ° r« 


y + -\jy 2 +v4 2 etc' 


where 


/l 2 = (x-;c') 2 +r 2 

Now for the first integral in equation (A7) 

4cv-, o )W(^-^) 2+ Cf-, o ) 2+ -^ 

471 e n L J 


+ _y 0 )W * 2 +(y-j 0 ) 2 + ^ 2 


47te 


471 €f 


Vv 2 + bv -z 2 + — tn\ 


2v/v 2 + t>v — z 2 + 2v + 6 


- zsm 


-l 


v* 


bv-2z 


IvlVt’ 2 + 4z 2 


where 6 = 2 (y 0 - _>-) and _>•’ < y 0 always 


V A = G" - A'o )-•- a/ (^0 - ^) 2 +z2 +CF-70) 2 

v « =0 ; -3 ; o) + V jr2 +r2 +G'-y 0 ) 2 

where we have used 380.31 1, 261.01, 380.001, and 380.1 1 1 from reference 87. The second term of 
equation (A7) is the negative of equation (A9) with yo = 0. 

On the upper surface we have 


r 9 - xa x + y f a v + z 0 a z 

r =xd x +yd y +zd,\r-r'\ = [(x-x'f +{y-y'f + {z-z 0 f] /2 

ds' = dx'dy' 


(A6) 

(A7) 

(A8) 

(A9) 


(A 10) 

(All) 
(A 12) 
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then 


«,(')= fA 

4 71 6n 


rr 


dy'dx' 


/2 


(A13) 


By observation of equation (A6), we see equation (A 1 3) only differs by a negative sign and the 
replacement of z 2 with (z-z 0 ) 2 : thus only minor modifications are necessary on equation (A9), to evaluate 
equation (A 13). The superposition of the results of equation (A9) and (A 13) constitute the complete 
potential of the two sheets of charge. The valid range is y 0 <y , and all x and z. 
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